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Abstract 

We have calculated the mass distribution (da/dMx) as observed by UA8 
Collaboration in the inclusive reaction pp pXp at y/s = 630 GeV, using 
the Interacting Gluon Model (IGM) with Double Pomeron Exchange (DPE) 
included. The only new parameter is the IP-IP cross section, which we can 
extract from fitting experimental data. We compare our results with the val- 
ues obtained in the UA8 study. Assuming a constant Pomeron-Pomeron total 
cross section (djpjp = Imb), we make predictions for da/dMx at Tevatron 
and LHC energies. 



I. INTRODUCTION 

After ten years of work at HERA, an impressive amount of knowledge about the Pomeron 
has been accumulated, especially about its partonic composition and parton distribution 
functions. Less known are its interaction properties. Whereas the Pomeron-nucleon cross 
section has been often discussed in the literature, the recently published data by the UA8 
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Collaboration [1] have shed some light on the Pomeron-Pomeron interaction. In [1] the 
Double Pomeron Exchange (DPE) cross section was written as the product of two flux 
factors with the IP-IP cross section, (Jjpp^ being thus directly proportional to this quantity. 
This simple formula relies on the validity of the Triple-Regge model, on the universality of 
the Pomeron flux factor and on the existence of a factorization formula for DPE processes. 
However, for these processes the factorization hypothesis has not been proven and is still 
matter of debate [2-5]. In [6] it was shown that factorizing and non- factorizing DPE models 
may be experimentally distinguished in the case of dijet production. 

Fitting the measured mass spectra allowed for the determination of ajpp and its depen- 
dence on Mx-i the mass of the diffr active system. The first observation of the UA8 analysis 
was that the measured diffractive mass (Mx) spectra show an excess at low values that can 
hardly be explained with a constant (i.e., independent of Mx) CFipip. Even after introducing 
some mass dependence in ctjpjp they were not able to fit the spectra in a satisfactory way. 
Their conclusion was that the low Mx excess may have some physical origin like, for example 
glueball formation. 

Although the analysis performed in [1] is standard, it is nevertheless useful to confront 
it with other, also successful, descriptions of the diffractive interaction. One of them is the 
one provided by the Interacting Gluon Model (IGM) [7] . This model describes only certain 
aspects of hadronic collisions, related to energy fiow and energy deposition in the central 
rapidity region. It should not be regarded as an alternative to a field-theoretical approach to 
diffractive amplitudes, but rather as an extension of the naive parton model. The reason for 
using it here is that it may be good enough to account for energy fiow in an economic way. 
The deeper or more subtle aspects of the underlying field theory probably (this is our belief) 
do not manifest themselves in energy flow, but rather in other quantities like the total cross 
section. Inspite of its simplicity, this model can teach us a few things and predict another 
few. This is encouraging because in the near future new data about DPE from CDF will be 
available [8]. In this work we would like to address the UA8 data with the IGM. As it will 
be seen, according to our analysis the low mass behavior of the Mx spectra can be explained 
with a constant IP-IP cross section. 

II. THE IGM 

The IGM has been described at length, especially in [7] and more recently in [11]. In the 
past we have successfully modified the IGM in such way as to include in it hadronic single 
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diffractivc dissociation processes [9-12] and applied it to hadronic collisions and HERA- 
photoproduction data [M]^ distributions and leading particles spectra of p, n, K and J/ip). 

The main idea of the model is that nucleon-nucleon collisions at high energies can be 
treated as an incoherent sum of multiple gluon-gluon collisions, the valence quarks playing a 
secondary role in particle production. While this idea is well accepted for large momentum 
transfer between the colliding partons, being on the basis of some models of minijet and 
jet production (for example HIJING [13]), in the IGM its validity is extended down to low 
momentum transfers, only slightly larger than Aqcd- At first sight this is not justified 
because at lower scales there are no independent gluons, but rather a highly correlated 
configuration of color fields. There are, however, some indications coming from lattice QCD 
calculations, that these soft gluon modes are not so strongly correlated. One of them is the 
result obtained in [14], namely that the typical correlation length of the soft gluon fields is 
close to 0.3 fm. Since this length is still much smaller than the typical hadron size, the gluon 
fields can, in a first approximation, be treated as uncorrelated. Another independent result 
concerns the determination of the typical instanton size in the QCD vaccum, which turns 
out to be of the order of 0.3 /m [15]. As it is well known (and has been recently applied to 
high energy nucleon-nucleon and nucleus-nucleus collisions) instantons are very important 
as mediators of soft gluon interactions [16]. The small size of the active instantons lead to 
short distance interactions between soft gluons, which can be treated as independent. 

These two results taken together suggest that a coUision between the two gluon clouds 
(surrounding the valence quarks) may be viewed as a sum of independent binary gluon-gluon 
collisions, which is the basic idea of our model. Developing the picture above with standard 
techniques and enforcing energy-momentum conservation, the IGM becomes the ideal tool to 
study energy flow in high energy hadronic collisions, in particular leading particle production 
and energy deposition. Confronting this simple model with several and different data sets 
we obtained surprisingly good agreement with experiment. 

As indicated in the recent literature [2-6], one of the crucial issues in diffractivc physics 
is the possible breakdown of factorization. As stated in [3] one may have Regge and hard 
factorization. Our model does not rely on any of them. In the language used in [3] , we need 
and use a "diffractivc parton distribution" and we do not really need to talk about "flux 
factor" or "distribution of partons in the Pomeron" . Therefore there is no Regge factorization 
implied. However, we will do this connection in Eq. (13), in order to make contact with the 
Pomeron pdf's parametrized by the HI and ZEUS collaborations. As for hard factorization, 
it is valid as long as the scale /i is large. In the IGM, as it will be seen, the scale is given by 



3 



/i^ = xys^ a number which sometimes is larger than 3 — 4 GeV^ but sometimes is smaller, 



part of the phase space we are inside the validity domain of hard factorization, but very 
often we are outside this domain. From the practical point of view, Eq. (7), being defined 
at a semihard scale, relies on hard factorization for the elementary gg gg interaction, uses 
parton distribution function extracted from DIS and an elementary cross section CTgg taken 
from standard pQCD calculations. The validity of the factorizing-like formula Eq. (6) of our 
paper is an assumption of the model. In fact, the relevant scale there is mg — ^^qcd 
strictly speaking, there are no rigorously defined parton distributions, neither elementary 
cross sections. However, using Eq. (6) has non-trivial consequences which were in the past 
years supported by an extensive comparison with experimental data. 



Double Pomeron exchange processes, inspite of their small cross sections, appear to be 
an excellent testing ground for the IGM because they are inclusive measurements and do 
not involve particle identification, dealing only with energy fiow. In what follows, we briefiy 
mention our main formulas. For further discussion we refer to the works [9-11]. 



In Fig. 1 we show schematically the IGM picture of a double Pomeron exchange event 
in a proton- antiproton coUision. The interaction follows then the usual IGM [7] picture, 
namely: the valence quarks fly through essentially undisturbed whereas the gluonic clouds 

of both projectiles interact strongly with each other (by gluonic clouds we understand a sort 
of "effective gluons" which include also their fluctuations seen as qq sea pairs). The proton 
(antiproton) looses fraction x {y) of its original momentum and gets excited forming what 
we call a leading jet carrying Xp = 1 — x {xp = 1 — y) fraction of the initial momentum. 

According to the IGM [1 1] , the probability to form a flreball carrying momentum fractions 
x and y of two colliding hadrons (see Fig. 1) is given by: 



going down to values only slightly above Ag^^. When the scale is large (/i^ > Pt^^^) we 
employ Eq. (7) and when it is smaller (ttIq < fi^ < Pr^-^) we use Eq. (6). Therefore, in 



III. DOUBLE POMERON EXCHANGE 



A. Kinematics 
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^ {y'){x-{x)r + {x'){y-{y)r-2{xy){x-{x)){y-{y)) I (1) 



X exp 
where 

D,y^{x'){y')-{xyy ; {x^y^ = dx' x'^ dy' y'^u;{x',y'), (2) 

JO Jo 

with xo being a normahzation factor defined by the condition that 

[ dx [ dyx{x,y)0{xy-Kl^iJ^l (3) 

JO JO 

with Kmin = ^ being the minimal inelasticity defined by the mass mo of the lightest possible 
central gluonic cluster. 

As it can be seen from (1), the probability that the incoming hadrons release an energy 
of ^/xys is a two-dimensional gaussian function of x and y with maximum governed by the 
momenta (x) and {y) (particular cases of (2)), which, in turn, depend on the integration 
limits Xmax s-nd ymax- By reducing these maximal value we select events in which the energy 
released by the proton and by the antiproton is small (i.e., Mx is small) and at the same time 
two rapidity gaps will be formed. This is how we define our "kinematical Pomeron" : a set 
of gluons belonging to the proton (or antiproton) carrying altogether a small fraction of the 
parent hadron momentum. The function u!{x', y') will be discussed below. In the formulation 
of the IGM Xmax — Umax = 1 if Only non-diffractive processes are present. In the model, 
diffraction (double Pomeron exchange) means reducing ymax {xmax a-nd ymax)- Although this 
procedure is somewhat arbitrary and we could choose any small number for the integration 
hmits, this freedom of choice is dramatically reduced if we use single diffractive events (SPE) 
as a guide. In [9] we have shown that the choice leading to the best description of diffractive 
mass spectra is ymax — U- Actually, using this cut in (2) and some simple approximations 
(described in [9] and also in the appendix) we could obtain the analytical formula for the 
single diffractive mass spectrum: 



dNspE Xo 



dMl TTv/^Miji^j^^^p/^'^'n ^^^^ 



1 - cln 



M 



(4) 



where c is a constant (discussed below). This spectrum shape is in very good agreement 
with a wide body of data. Based on our previous success we shall assume here that in double 
Pomeron exchange we have Xmax = x and ymax = V and consequently Xmax Umax ^ xy = 

In [8] it has been conjectured that the ratio of two-gap to one-gap rates could be used to 
test QCD aspects of gap formation. We therefore calculate in the appendix, with the same 
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approximations, the analytical expression for the DPE mass spectrum, which turns out to 
be: 

2 



DPE 



Xo 



In 



s 

All 



exp < 



2 • 



1 - c'ln 



Ml 



d In 



a4 



(5) 



Both expressions (4) and (5) are dominated by the 1/M|^ factor. This garantees a 
pronounced fall with increasing diffractive masses, which is confirmed by the numerical 
calculations (which have no approximation) presented below. 



B. Dynamics 

The spectral function, uj{x',y'), contains all the dynamical input of the IGM. Their soft 
and semihard components are given by (cf. [17]): 

^'{x', y') = G{x') G{y') e {x'y' - K^^) , (6) 



y') = G{x') G{y') e (x'y' - ^] , (7) 

a{s) \ s J 

where G"s denote the effective number of gluons from the corresponding projectiles (approx- 
imated by the respective gluonic structure functions), 6"^^ and are the soft and semihard 
gluonic cross sections, PTmin is the minimum transverse momentum for minijet production 
and (7 = (Tjpip denotes the Pomeron-Pomeron cross section. 

In order to be more precise, the function G{x') {G{y')) represents the momentum dis- 
tribution of the gluons belonging to the proton (antiproton) subset called Pomeron and x' 
{y') is the momentum fraction of the proton {antiproton) carried by one of these gluons. We 
shall therefore use the notation G{x') — G]p{x'). This function should not be confused with 
the momentum distribution of the gluons inside the Pomeron, fg/ip{P). 

The Pomeron for us is just a collection of gluons which belong to the diffracted proton 
(antiproton). In our previous works we have assumed that these gluons behave like all other 
ordinary gluons in the proton and have therefore the same momentum distribution. The 
only difference is the momentum sum rule, which for the gluons in IP is: 

['dx'x'Gjp{x') = pd (8) 
Jo 
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where pd = 0.05 (see [9,12]) instead of p = 0.5, which holds for the entire gluon population 
in the proton. 

In order to make contact with the analysis performed by HERA experimental groups we 
consider two possible momentum distributions for the gluons inside IP. A hard one: 

fi/AP) = ah{l- (3) (9) 
and a "super-hard" (as it is called in [5]) or "leading gluon" (as it is called in [18]) one: 

rgfAP) = ashP'{l - (3)'-' (10) 

where f3 is the momentum fraction of the Pomeron carried by the gluons and the superscripts 
h and sh denote hard and superhard respectively. The constants and ash will be fixed by 
the sum rule (8). In the past [10], following the same formalism, we have also considered a 
soft gluon distribution for the Pomeron of the type 

(1 - 0Y 

/.X/5)=6^-^ (11) 

but we found that this "soft Pomeron" distribution was incompatible with the single diffrac- 
tive mass spectra measured at HERA [19]. This Pomeron profile was also ruled out by other 
types of observables, as concluded in Refs. [20] and [21]. 

We shall use the Donnachie-Landshoff Pomeron flux factor, which, after the integration 
in the t variable, is approximately given by [5]: 

fjp/p{xjp) ~ Cx]^'^" ~ C — (12) 

X]p 

where X]p is the fraction of the proton momentum carried by the Pomeron and the normal- 
ization constant C will be fixed later, also with the help of (8). Noticing that P — the 
distribution G]p{y) needed in eqs. (6) and (7) is then given by the convolution: 

G'/'iy) - fjP/M f^;^^) (13) 

Jy X]p X]p 

We shall use also the "diffractive gluon distribution" given by: 

(1 - v)^ 

GAv) = « ^ — (14) 

y 

where a is fixed by the sum rule. With (14) we could obtain a very good description of 
diffractive mass spectra [9,10]. Therefore we shall keep using it here. 
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We are implicitly assuming that all gluons from p and p participating in the collision 
(i.e., those emitted from the upper and lower vertex in Fig. 1) have to form a color singlet. 
Only then two large rapidity gaps will form separating the diffracted proton, the Mx system 
and the diffracted antiproton, which is the experimental requirement defining a DPE event. 

We can now calculate the diffractive mass distribution Mx using the x(a;, y) function by 
simply performing a change of variables: 

1 da 



adMx 



J^dx J^dy x{x, y) S [Mx - ^/xys\ 9 (xys - mg) 

ml) 



2M 



^2 dx \ 



x,y 



Ml 



9{Ml 



(15) 



IV. NUMERICAL RESULTS AND DISCUSSION 

We start evaluating Eq. (15) with the inputs that were already fixed by other applications 
of the IGM [9,10], namely, (14) with — 0.05. In Fig. 2 we show the numerical results 
for DPE mass distribution. We have fixed the parameter a (= (Tpp) appearing in Eq. (6) 
and Eq. (7), to 0.5 mfe (solid fines) and 1.0 m6 (dashed fines). In both cases our curves were 
normafized to the "AND" (Fig. 2a) and "OR" (Fig. 2b) data samples of [1] . We emphasize 
that, in this approach, since we have fixed all parameters using previous data on leading 
particle formation and single diffractive mass spectra, there are no free parameters here, 
except ajpjp. 

As one can see from the figures, in our model we obtain the fast increase of spectra in the 
low mass region without the use of a Afx dependent P-IP cross section and this quantity 
seems to be approximately apjp ^ 0.5 m6. This is the main message of this note. 

In order to investigate the sensitivity of Eq. (15) to the value of pd-, we show in Fig. 3 
with arbitrary units, diffractive mass spectra obtained with p^ = 0.025 (dash-dotted line), 
0.05 (solid line) and 0.1 (dashed line). The three curves have the same normalization and we 
observe that, as pd increases our mass spectrum becomes softer. It is interesting to remark 
that, in the actual calculations, the parameter pa appears always divided by ajpjp in the 
computation of the moments (2). In fact, they form one single parameter. Assuming that 
the Pomeron profile is universal, we could disentangle one from the other, fitting p^ from the 
analysis of previous data [9,10] and now extracting ctjpjp from the the UA8 data. 

We next replace (14) by the convolution (13) to see which of the previously considered 
Pomeron profiles, hard or superhard, gives the best fit of the UA8 data. In doing so, we shall 
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keep everything else the same, i.e., = 0.05 and ajpjp = 0.5 m6. In Fig. 4a we compute 
(13) and show yG]p{y) for Eq. (14) (sohd hne), the hard distribution (9) (dashed hue) 
and the super hard (10) (dash-dotted hne). For the sake of comparison. Fig. 4b shows the 
corresponding diffractive mass spectra normahzed to the unity with same notation. We see 
that, for harder Pomeron profiles we "dig a hole" in the low mass region of the spectrum. 

In Fig. 5, we repeat the fitting procedure used in Fig. 2 for the Pomeron profiles shown 
in Fig. 4. We fix pd = 0.05, ctjpjp = O.bmb. Solid, dashed and dash-dotted lines represent 
respectively Eq. (14), hard and superhard Pomerons. Note that the solid lines are the same 
as in Fig. 2. Looking at the figure, at first sight, we might be tempted to say that Eq. (14) 
gives the best agreement with data and a somewhat worse description can be obtained with 
the hard Pomeron (in dashed lines), the superhard being discarded. However, comparing 
the dashed lines in Fig. 2 and Fig. 5 and observing that they practically coincide with each 
other, we conclude that the same curve can be obtained either with (14) and a^jp — 1.0 mb 
(dashed hne in Fig. 2) or with (13), (9) and app — 0.5 mb (dashed line in Fig. 5). In other 
words we can trade the "hardness" of the Pomeron with its interaction cross section. The 
following two objects give an equally good description of data: i) a Pomeron composed by 
more and softer gluons and with a larger cross section and ii) a Pomeron made by fewer, 
harder gluons with a smaller interaction cross section. We have checked that this reasoning 
can be extended to the superhard Pomeron. Although, apparently disfavoured by Fig. 5 
(dash-dotted fines), it might still fit the data provided that a pip < 0.25 mb. Given the 
uncertainties in the data and the limitations of the model, we will not try for the moment 
to refine this analysis. It seems possible to describe data in a number of different ways. We 
conclude then that nothing exotic has been observed and also that the Pomeron- Pomeron 
cross section is bounded to be ajpjp < l.Omb. 

In Fig. 6a we compare our predictions for da/dMx {mb/GeV) at Tevatron {^/s = 2TeV) 
and in Fig. 6b for LHC {^/s = lATeV) assuming an Mx -independent ajpjp = l.Omb (and 
using (14)) with predictions made by Brandt et al. [1] for two values of effective Pomeron 
intercepts («(0) = 1 + e), e = 0.0 and 0.035. 

Although the normalization of our curves is arbitrary, the comparison of the shapes 
reveals a striking difference between the two predictions. Whereas the points (from [1]) show 
spectra broadening with the c.m.s. energy, we predict (solid lines) the opposite behavior: as 
the energy increases we observe a (modest) narrowing for da/dMx- This small effect means 
that the diffractive mass becomes a smaller fraction of the available energy ^/s. In other 
words, the "diffractive inelasticity" decreases with energy and consequently the "diffracted 
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leading particles" follow a harder xp spectrum. Physically, in the context of the IGM, this 
means that the deposited energy is increasing with ^/s (due the minijets) but it will be 
mostly released outside the phase space region that we are selecting. 

We are not able to make precise statements about the diffractive cross section (in par- 
ticular about its normalization) with our simple model. Nevertheless, the narrowing of 
daupE/dMx suggests a slower increase (with ^/s) of the integrated distribution aopE- We 
found this same effect [9] also for (Tspe- This trend is welcome and is one of the possible 
mechanisms responsible for the suppression of diffractive cross sections at higher energies 
relative to Regge theory predictions. 

In Fig. 7 we show the ratio R{Mx) defined by: 

1 danpE 

RiMx) = Z:. (16) 

CSPE dMx 

This quantity involves only distributions previously normalized to unity and does not di- 
rectly compare the cross sections (which are numerically very different for DPE and single 
diffraction). In R the dominant factors cancel, as suggested by the comparison be- 

tween (4) and (5), and we can better analyse the details of the distributions which may 
contain interesting dynamical information. The most prominent feature of Fig. 7 is the 
rise of the ratio with Mx, almost by one order of magnitude in the mass range considered. 
This can be qualitatively attributed to the fact that, in single diffractive events the object 
X has larger rapidities than the corresponding cluster formed in DPE events. As a conse- 
quence, when energy is released from the incoming particles in a SPE event, it goes more 
to kinetic energy of the X system (i.e., larger momentum Px and rapidity Yx) and less 
to its mass. In DPE, although less energy is released, it goes predominantly to the mass 
Mx of the difractive cluster, which is then at lower values of Yx- In order to illustrate this 
behavior, we show in Fig. 8 the rapidity distributions of the X system (which has Mx)- 
All curves are normalized to unity and with them we just want to draw attention to the 
dramatically different positions of the maxima of these distributions. The solid and dashed 
hnes show 1/ada/dYx for DPE (curves on the left) and SPE (curves on the right) computed 
at ^/s — 630 GeV and a/s — 2000 GeF, respectively. We can clearly observe that DPE and 
SPE rapidity distributions are separated by three units of rapidity and this difference stays 
nearly constant as the c.m.s. energy increases. The location of maxima in 1/ada/dYx and 
their energy dependence are predictions of our model. 

To summarize: we have further developed our model for hadronic collisions and included 
double Pomeron exchange events. With only one new parameter, ajpjp, we could fit the data 
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recently published by the UA8 collaboration and make predictions for the DPE mass spectra 
at Tevatron and LHC energies. Our main conclusion is that ajpjp ~ 0.5 mb and constant 
with Mx is favored by experimental data. We predict that the ratio between (normalized) 
double ]P-]P exchange and single diffractive mass distributions grows with Mx- 



V. APPENDIX 



In what follows we shall often make use of our kinematical constraint between Xm{= x^ax) 

and ym{= Umax)- 

Mx 

XruVm =Xy = (17) 

s 

The moments of u!^{x' , y') and u!^{x', y') are given by: 

{x^y^s = fjdx'x'^ fjdy'y'^u'{x\y') = j^.dx' x'^ dy' y'^ u\x\y') (18) 



{x''ynH= / dx'x"' / dy'y"^uj''{x',y')^ L4 dx' x'^ ^ dy' y"^ uj"" {x\y') (19) 

M ^ s 

X 

In order to obtain analytical expressions we shall, in the following set uj^{x',y') = 
because at the relevant energies hard processes are not yet dominant. We shall keep only 
the low X dominating factor of the gluon distribution: 

Gipix') = i (20) 
x' 

We shall neglect the x' and y' dependence of the cross sections and assume that: 

= c (21) 
With all these approximations (18) can be rewritten as: 



{x'^yn 



= dx' x"" dy' 2/"" u{x', y') 9 I x'y' - ^ j 

= c dx'x'--^ dy'y'-^ 9 [x'y' - ^) 

^dx'x'^-' dy'y'^-\ (22) 
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By symmetry we have: 
(y) ~ c- 
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B. Case n=2, m=0 

dy' J / / , ( VmS 



\x ) = c dx X L2 —- = cl2 dx X m 

- ^ SV-m \ ^ 



s 

Again, by symmetry, we have: 
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, 2, 1 /Mix 1 , , 1 



C. Case n=l, m=l 



(xy) 



rym ^ 

dx dy 



1^ dx' 
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sym 


sx' 



S 



1^0 1 

— c • m 



. ^0 . 



c • 



'Ml- 



~ 



Inserting the approximate expressions for the moments into (1) we obtain: 

Xo _____ f {y^){x - {x) f + {x^){y - {y) f\ 



x{x,y) ^ 



{^'){y') 



exp 



2D 



xy 



or 
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x{x,y) 



Xo 



■ exp 



{x-{x)f {y-{y)f 



2(x2) 



The diffractive mass distribution will be given by: 
dN 



2{y') 



(30) 



J^dx J^dy x{x, y) 8 (m^ - ^ys) 6 \^xy 
1 dx 

- / — : 



.1 ^ \ xs ^ 



(31) 



In order to proceed further let us first notice that x^ and y^ in the formulas for < > 
above will have the meaning of the x and y because of the 6{M]^—xys) constraint. Therefore: 



(x) — a 



a 

y 



Xr, 



a 

X 



(32) 



la /M|\ l_h_ 
2y' \ s ' y y^ 



where 



a — c 



'Ml' 



In 



X 



la /M|~ 
2x ' 



X x^ 



It then follows that: 



exp 



exp 



- jx))^ _ {y-{y)f 
2{x^) 2(y2) 



2^ 



2^ 





{xy - a) 2 


— exp 


b 



Using once again that xy = 
{xy - af 

exp 



we arrive at: 

s 



— exp < 



exp < 



s \ s \ mf. 



1 - c[nl!^ 



cln ( 



The total x(x, y) is then: 



(33) 



(34) 



(35) 



(36) 
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leading to 



dN 



In 



Mi 



Xo 1 



m; 



exp 



1 — cln 



Mi 



2 ^ 



cln 



which is exactly Eq. (5). 
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FIGURES 




FIG. 1. IGM picture for a double Pomeron exchange process. 
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FIG. 2. IGM DPE diffr active mass distributions: Solid and dashed lines show the numerical 
results with ctjpjp equal to 0.5mb and 1.0 m6, respectively. Our curves were normalized to the 
"AND" (a) and "OR" (b) data samples of [1]. 
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FIG. 3. IGM DPE diffractive mass distributions, in arbitrary units, with pd = 0.025 
(dash-dotted line), 0.05 (sohd line) and 0.1 (dashed line). In all cases ajpjp = 0.5mb. 
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FIG. 4. a) Diffractive gluon distributions: solid line was calculated with Eq. (14), dashed and 
dash-dotted lines, calculated by Eq. (13), represent the "hard" and "super-hard" Pomeron profiles, 
given respectively by Eq. (9) and (10); b) Diffractive mass distributions, normalized to the unity, 
for the same cases showed in a). In all cases pd = 0.05 and ajpjp = 0.5mb. 
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FIG. 5. IGM DPE diffractive mass distributions: Solid, dashed and dash-dotted hnes fohow 



the same notation of Fig. 4. Our curves were obtained with pd = 0.05 and crjpjp = O.bmb and 
normahzed to the data samples of [1]. 
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FIG. 6. IGM predictions for da/dMx at Tevatron and at LHC with ctjpjp = l.Omft. Cross(-I-) 
and Cross(x) are predictions made by Brandt et al. [1] for two values of effective Pomeron intercepts 
a(0) = 1 + e. 
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FIG. 7. Ratio double/single Pomeron exchange mass distributions as a function of Mx- In both 
cases we have assumed ajpjp = 1.0 mb (for DPE processes) and apip = 1.0 mb (for SPE processes). 
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FIG. 8. Double and single Pomeron exchange normalized rapidity (Yx) distributions. In both 
cases we have assumed (Tjpjp = 1.0 mb (for DPE processes) and iTpjp = 1.0 mb (for SPE processes). 
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